First-principles study of diffusion behaviors of Cs and I in Cr coating* 
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Cs and I can migrate through the fuel-cladding interface and accelerate the cladding corrosion 
process induced by the Fuel Cladding Chemical Interaction. Cr coating is an important candidate in 
mitigating this chemical interaction. First-principles calculations were employed to investigate the 
diffusion behaviors of Cs and I in Cr bulk and Cr grain boundary, aiming to reveal the microscopic 
mechanisms for mitigating the interaction at fuel-cladding interface. The interactions between these 
two fission products and Cr coating were systematically studied, and the temperature-dependent 
diffusion coefficients of Cs and I in Cr were obtained using Bocquet’s oversized solute atoms model 
and Le Claire’s nine-frequency model, respectively. The results show that the migration barriers of 
Cs and I are significantly lower compared to that of Cr, and the diffusion coefficients of Cs and I are 
both more than 3 orders of magnitude larger than Cr self-diffusion coefficient within the temperature 
range of Generation IV fast reactors (below 1000 K), which shows the strong penetration ability of 
Cs and I. Meanwhile, Cs and I are more likely to diffuse along the grain boundary because of the 
generally low migration barriers, indicating that grain boundary serves as a fast diffusion channel 
for Cs and I. 
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I. INTRODUCTION 


Fuel Cladding Chemical Interaction (FCCI) is consid- 
ered as one of the major factors limiting the lifetime of 
fuel pins in fast reactors, especially for oxide fuel pins 
with stainless steel cladding [1], as severe oxidative corro- 
sions can occur on the interface of the fuel and the inner 
wall of the cladding. Some fission products, migrating 
through the fuel-cladding interface, can accelerate the 
FCCI-induced cladding corrosion process, and the main 
fission products involved are volatile Cs and I [1, 2] con- 
sidering their production yield and penetrating depth. 
As the cladding attack rate depends on diffusion of reac- 
tants at cladding surface, the diffusion of Cs and I in fuel 
pins is a key factor in FCClI-induced corrosion problem. 

Due to the complex and multiscale coupled physi- 
cal and chemical mechanisms of FCCI, current research 
mainly focuses on exploring materials that can mitigate 
FCCI through experimental studies, and a coating of 
buffer-getter materials on the cladding inner surface is 
considered as an effective method [1]. The buffer-getter 
materials including V, Nb, Cr, Zr, U, Ti and certain 
rare earth metals can reduce oxidation of the cladding 
due to their high affinity for oxygen [1]. Specifically, 
Cr coating features a high melting point, superior oxida- 
tive corrosion resistance [3, 4], good crack resistance [5- 
7], and mature preparation processes [8-10], making it 
an important candidate in mitigating FCCI, but the 
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underlying micro-mechanisms remain unknown. First- 
principles density functional theory (DFT) is recognized 
as a powerful tool to reveal the micro-mechanisms, espe- 
cially from the perspective of elemental diffusion behav- 
iors when coupled with harmonic transition state the- 
ory [11], as proved by numerous previous works for im- 
purity diffusion including its consequent effects in Ni [11— 
16], Fe [17-24], Al [25-27], SiC [28], etc. Recently, 
Yang et al. have carried out first-principles calculations 
of lanthanides diffusion in Cr and a-Fe to identify the 
FCCI-mitigating mechanisms for metallic fuel and fer- 
ritic/martensitic steel cladding system [22]. According 
to our knowledge, there are still few investigations of 
FCCl-involved elements for oxide fuels, in which Cs and 
I dominate. Besides, further investigation is warranted 
to comprehensively understand the effects of Cr coat- 
ing on the cladding inner surface functions in mitigating 
FCCI for oxides fuel pin systems. 

In this work, first-principles method was employed to 
investigate the diffusion behaviors of important fission 
products (Cs and I) in BCC Cr coating. Cs and I are 
so large in the Cr matrix that the diffusion mechanism 
by interstitial atoms [29-33] was deemed unrealistic. In- 
stead, vacancy-mediated mechanism is the main mode 
of diffusion for impurity atoms of large atomic size [34], 
such as Cs and I. Section 2 lays out the vacancy-mediated 
diffusion models and details of the first-principles calcu- 
lations. ‘The methodology to determine the inputs for 
diffusion models is also presented. Section 3 discusses 
the interactions between impurity atoms and the sur- 
rounding Cr atoms or vacancies. Diffusion coefficients 
of Cs and I in BCC Cr were obtained by applying the 
diffusion model, and possible diffusion paths of Cs and I 
at grain boundary (GB) were predicted based on results 
of migration barriers between different sites. 
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Fig. 1. (Color online) Illustration of diffusion models: (a) Le Claire’s model [35-37]; (b) OSA model [21]. 
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A. Diffusion models for vacancy-mediated solutes 


In this section, two diffusion models for determina- 
ss tion of impurity diffusion coefficients were introduced, 
namely Le Claire’s nine-frequency model [35-37] (re- 
ferred to as Le Claire’s model below) and the oversized 
solute-atom (OSA) model [21]. As shown in Fig. 1, im- 
purity atoms and matrix atoms are represented by or- 
ange and blue spheres, respectively, and vacancies by 
gray squares. ‘The numbers on the atoms and vacancies 
indicate their positions as the nearest neighbors (nn) to 
the impurity atom. The corresponding jump frequen- 
cies for different atom-vacancy exchange processes are 
also labeled in Fig. 1. The main difference between the 
two diffusion models lies in the configuration formed by 
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7s the impurity atom and its Inn vacancy, which will be 
zə described in details in the following. 

According to the Le Claire’s model [35-37], the self- 
diffusion coefficient (Dseif—aif¢) and the solute diffusion 
coefficient (Dsolute) in BCC lattice can be obtained by: 
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with a the lattice parameter, Æ? the vacancy-solute bind- 
ing energy when the vacancy is located on the Inn site 
of the solute, kg the Boltzmann constant and T the 
temperature in K. fo is the correlation factor for self- 
diffusion, and is a constant equal to 0.727 for diffusion 
in BCC crystals. The correlation factor for solute diffu- 
sion f, can be expressed as [35-37]: 
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where F = 0.512. This expression applies to the case 
where the impurity atom is close in size to the ma- 
trix atom, with approximation wg = wo. Based on the 
transition-state theory [35], the average atomic jump fre- 
quency w; can be written with respect to the migration 
energy (£7) as: 
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in which v; and £7” are the attempt frequency and cor- 
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responding migration barrier for jump i, respectively. vi 
, _ m3N-3 7 /m3N-4 T 
is defined as v; = | [j= 1/1} 4 | 
normal vibration frequencies of the initial state (I) of 
migration, and v? is the nonimaginary normal frequen- 
cies of the transition state (T). The vacancy concentra- 


tion at equilibrium at temperature T is given by Cy = 


in which vi is 
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exp(—7—p )exp(;2), with the vacancy formation energy 
19 Ef = E! —(N—1)E®/N and the vacancy formation en- 
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Fig. 2. (Color online) Illustration of simulation models for Cr matrix and GB. (a) 4 x 4 x 4 supercell containing 128 atoms 
representing Cr matrix; (b) Cr ©5(210) GB with 76 atoms, where the atomic sites at GB interface are labeled by the numbers 


(0 for interstitial site, 1 & 2 for substitutional sites). 


which EŻ and EË are the corresponding energy of the 
initial state structure and the bulk structure containing 
N — 1 and N atoms, respectively, and vË is normal vi- 
bration frequencies of the bulk structure (B). 

The Le Claire’s model is the commonly used diffusion 
model which describes well the diffusion of substitutional 
solutes. Nevertheless, due to the strong attraction of the 
OSA to the vacancy (as illustrated in Fig. 1 (b)), the dif- 
fusion of OSA cannot be described by Le Claire’s model. 
Bocquet et al. proposed a diffusion model applicable to 
OSA, i.e. the OSA model, and the expression for OSA 
diffusion coefficient is given by [21]: 


2 


Dosa = igi maf, 


(5) 
where A is the distance between the vacancy and its 
nearest atom, i.e. A = V3a/2 in BCC crystals, rmy 
is the frequency attached to a macrojump (MJ), and f 
is the correlation factor. A macrojump is formed by 2 
elementary displacements: the OSA located on a substi- 
tutional site (S) is pushed onto a transitional site (T), 
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The hypothesis of the “one-shot” approximation is that, 
after the dissociation of the OSA and its nearby vacancy, 
only one jump is needed for the vacancy to return close 
to the OSA. It is proved available when the vacancy- 
solute interaction is very strong at lnn, which is exactly 
the case of OSA diffusion. 


Our preliminary calculation show that Cs combines 
with its Inn vacancy and forms a configuration of “V /2 
159 + Cs + V/2” (V for vacancy), corresponding to Fig. 1(b). 
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130 i.e., S — T, and then moved onto a S site from T site 
i31 (T —> S). In this way the frequency attached to a MJ is 
i32 defined as: [yy = Usrlrs/(Usr F Trs) with Por = 
b b b 
133 8Cy [Bwreap( <7) + Bw exp( zr) + wan exp(q—r)] and 
134 Irs = x (3w3 + 3w3 + wg ). EP is the binding energy 
135 between the solute and its i*” nn vacancy (positive value 
136 for attraction, and vice versa), and can be calculated by: 


(6) 


i3s The terms on the right side of the equation represent the 
energies of different configurations, and the numbers of 
Cr atoms, vacancies, and impurity atoms in each configu- 
ration are (127, 1, 0), (127, 0, 1), (126, 1, 1), and (128, 0, 
0), respectively. The correlation factor f for OSA diffus- 
ing in BCC crystals equals to 1+ Q?°°, in which Q?°° 
is the average cosine between a T — S jump vector and 
145 the next S > T jump in the BCC lattice. Bocquet et 
146 al. obtained the exact formula of Q?°° by a double 
147 Laplace and Fourier transform of the transport equation 
148 for the vacancy. But for simplification, here we applied 
149 the formula with the “one-shot” approximation [21]: 
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ico In contrast, I does not exhibit strong attraction to Inn 
161 V, corresponding to Fig. 1(a). According to the applying 
162 conditions for different diffusion models, the OSA model 
163 is used for the calculation of Cs diffusion coefficients, and 
164 the Le Claire model is used for I in this study. 
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B.  First-principles calculation 


To investigate the diffusion behaviors of Cs and I in 
BCC Cr coating, the density functional theory (DFT) 
calculations are performed with the Vienna ab initio sim- 
ulation package (VASP) [39-41] using the projector aug- 
mented wave (PAW) [42] pseudopotentials and Perdew- 
Burke-Ernzerhof (PBE) [43] parametrization. Razu- 
movskiy et al. found that disordered local-moment cal- 
culations for pure Cr in the paramagnetic state will re- 
duce to a nonmagnetic solution [44]. A 4 x 4 x 4 super- 
cell containing 128 atoms was adopted for the simulation 
of Cr matrix, and the ©5(210) GB supercell consisting 
of 76 sites for simulation at Cr GB, as shown in Fig. 2. 
The Brillouin area was sampled with the 4 x 4 x 4 and 
5 x 5x 1 k-point meshes for Cr matrix and GB, respec- 
tively. The total energy and force convergence criteria 
were set to 1.0 x 107° eV and 0.01 eV/A, respectively, 
with a cutoff energy of 400 eV. All transition states and 
migration barriers were obtained by applying the climb- 
ing image Nudged Elastic Band (CNEB) method, using 
3 intermediate images. To determine the diffusion co- 
efficients in Cr matrix, phonon calculations were also 
performed considering only the vibrational modes of the 
migrating atoms |22, 45-48]. For both the search of tran- 
sition states and phonon calculations, the force conver- 
gence criterion was reset to 0.001 eV/A, as the results 
are highly sensitive to this parameter. 


HI. RESULTS AND DISCUSSION 


A. 


Impurity diffusion in Cr matrix 
1. Atomic structure and electron distribution in matrix 


The interaction between atoms can change the forces 
applying on the impurity, thus is an important factor 
affecting the migration behaviors. Fig. 3 illustrates the 
change in bond lengths by a substitutional impurity. It 
can be seen that both Cs and I exert a repulsive effect 
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on the surrounding Cr atoms, because the Bader radii of 23; 
Cs and I shown in Table 1 are both larger than that of 23. 


Cr. The effect of Cs is more significant, with the bond 
lengths of its Inn and 2nn increasing both by more than 
0.1 A. 

To better understand the effect of impurities, the 
electronic structure with impurities was further ana- 
lyzed. Using the Bader decomposition method [49], 
Bader charge, volume and radius [50] were listed in Ta- 
ble 1. For the Cr atoms, it is found that the Bader charge 
and atomic volume of the 127 Cr atoms are almost the 
same, so their average values are given directly in the ta- 
ble. It is observed that both Cs and I gain electrons from 
the matrix. In addition, the non-metallic I atom gains 
reasonably more electrons than the metallic Cs atom. 
The atomic volumes of Cs and I are almost the same 
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when they occupy a substitutional site in Cr matrix, so 
it can be expected that there should be similarities in 
their migration behaviors. 
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Fig. 3. (Color online) Calculated bond lengths (A) in the Cr 
matrix for (a) pure matrix, (b) matrix with a substitutional 
Cs atom, (c) matrix with a substitutional I atom. 


TABLE 1. Calculated Bader charge,volume and radius of 
atoms in Cr matrix with Cs or I substitution. Cr(Cs) denotes 
the Cr atom in Cs-doped Cr matrix, and Cr(I) means that in 
I-doped one. 


Cs Cr(Cs) I Cr(I) 

Bader charge/e7 9.369 5.997 7.784 5.994 
Atomic volume/A® 18.778 11.346 19.056 11.344 
Atomic radius/A 1.649 1.394 1.657 1.394 


The binding effect between impurity atoms and their 
surrounding Cr atoms is weaker compared to that be- 
tween Cr atoms, which can be observed from the charge 
densities between the impurity atoms and their Inn Cr 
atoms shown in Fig. 4. It can be seen that the charge 
density around the impurity atoms decreases signifi- 
cantly compared to the corresponding profile for the pure 
Cr system (Fig. 4(a)), which is in agreement with the ef- 
fect of increasing bond length in Fig. 3. Besides, the 
presence of the impurity atoms also increases the charge 
density between the Cr atoms in the diagonal direction, 
indicating an enhanced interaction between them. Taken 
together, the Cr atoms near the impurity atoms tend to 
move away from the impurity atoms and have a weak 
binding interaction with the impurity atoms, favoring 
the diffusion of the impurity atoms in the matrix. 

In addition, the binding effect of Cr atoms to Cs is 
slightly weaker than that to I. Comparing the results of 
Cs and I (Fig. 4(b) & (c)), the charge density distribu- 
tions are similar, but in terms of equipotential lines, the 
charge density around Cs is more dilute and the interac- 
tion with Cr atoms is weaker than that of I. This is also 
consistent with the fact that the Cs-Cr bonds are longer 
than the I-Cr bonds in Fig. 3, and can also be observed in 
Fig. 5. There is no hybridization peak between Cs and Cr 
atoms, indicating the metallic bonds; several hybridiza- 
tion peaks appear between I and Cr atoms, indicating 
covalent-like type bonds. Specifically, the d electrons of 
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Fig. 4. (Color online) Calculated charge density distribution (electron/Bohr®) in the (101) plane for (a) pure Cr matrix, (b) 
Cr matrix with a substitutional Cs atom, (c) Cr matrix with a substitutional I atom. 


Cr and the p electrons of I are mainly involved in hy- 


Cr and I are mainly involved in hybridization between -2 
and 0 eV. Overall, Cs and I have similar effects on the Cr 
matrix, with slightly lower interactions around Cs than 


252 I. 


253 2. Diffusion process and diffusion coefficient in matrix 


In order to investigate the effect of impurity atoms on 
neighboring vacancies, the binding energies of impurity 
atoms to neighboring vacancies were calculated accord- 
ing to Eq. (6) and the results from 1nn to 4nn are shown 
in Fig. 6. The binding energy involved the 5nn vacancy 
is not presented, because the 5nn vacancy slides spon- 
taneously into the 1nn position, i.e., the Cr atom in the 
lnn position of the impurity slides spontaneously into 
the 5nn position. Both impurity atoms have an attrac- 
tive effect on the vacancies from 1nn to 4nn. The binding 
energy of Cs to the Inn vacancy is up to 2.46 eV, which 
is a possible reason for the formation of “V/2 + Cs + 
266 V/2” configuration. The relatively low binding energy 
between I and its Inn vacancy may be the reason why I 
remains stable at the substitution site. ‘The binding en- 
ergies of Cs atom to vacancies are generally higher than 
those of I atom, especially for the 1nn and 2nn vacancies. 
But for the 3nn and 4nn vacancy, the binding energies 
are almost the same, indicating that the impurity atoms 
mainly affect the Inn and 2nn vacancies. 

Migration energy and the attempt frequency required 
for the vacancies near the impurity to move onto differ- 
ent sites are also important factors affecting the difficulty 
of migration processes. ‘The results in Table 2 indicate 
278 that the impurity atoms are favored in the competition 
279 for vacancies with the surrounding matrix atoms. The 
230 Migration barrier of vacancies moving away from the im- 
231 purity atoms (£3", E}, Ep) are generally higher than 
2832 those in the direction close to the impurity atoms (E7’, 
233 Ep, Eg), and the energy required for the vacancies lo- 
234 cated in the 3nn of the impurity atoms to migrate to the 
285 Inn is lower than that required for self-diffusion in the 
286 pure Cr system (0.92 eV shown in Fig. 7). 
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bridization between -9 and -3 eV, and the d electrons of 2ss 


In terms of 309 


the attempted frequency, the attempt frequencies of va- 
cancies moving away from the impurity atoms (v3, V3, 
Vs) are also generally lower than the attempt frequencies 
near the impurity atoms (v4, v4, vg). This indicates that 
vacancies are more likely to be attracted to stay near the 
impurity atoms, which is also consistent with the results 
we obtained in Fig. 6, and this attractive effect to va- 
cancies facilitates the diffusion of impurity atoms. 
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TABLE 2. Vacancy migration energies and attempt frequen- 
cies for vacancies nearby the Cs or I atom. 


Cs I 

ES 7 0.36 
EY 2.70 1.66 
EY 1.19 0.96 

- 2.66 ie 

n 0.35 0.59 
p” J 1.09 
EX Í 0.59 
V2 / 2.85 
V3 5.06 1.22 
V4 12.80 10.12 
Va 9.79 10.04 
V4! 12.95 15.12 
V5 7 8.74 
U6 J 12.63 


Both Cs and I can diffuse easily in Cr matrix because 
of their low migration barriers. As illustrated in Fig. 7(a) 
& (b), Cs will spontaneously move from the substitu- 
tional site to the middle point of the two lattice sites, 
while I should overcome a migration barrier of 0.36 eV 
to move to another substitutional site. Although the 
migration of I is not spontaneous, the migration barrier 
is lower than that for Cr atoms (0.92 eV according to 
our calculations). As a result, the diffusion of Cs and I 
are both easier than the Cr self-diffusion. Nevertheless, 
305 the results of migration barriers do not necessarily imply 
306 that Cs diffusion is faster than I, as the displacement of 
307 I that occurs by diffusion is /3a/2, which is twice of the 
30s displacement of Cs. 
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Before determination of the diffusion coefficients of Cs 


(a) Pure Cr 


(b) Cr (Cs) 


Density of States (electron/eV) 


Energy (eV) 


Density of States (electron/eV) 


Energy (eV) 


Fig. 5. (Color online) Density of states (DOS) with respect to the energy: (a)&(d) Cr in pure Cr matrix, (b)&(e) Inn-Cr of 
a substitutional Cs or I atom, and (c)&(f) substitutional Cs and I atoms, respectively. Main hybridization peaks are marked 


with light green bars. 
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Fig. 6. (Color online) Calculated impurity-vacancy binding 
energies. 


and I in Cr, approximations should be made for the phe- 
nomenon of spontaneous sliding of the 5nn vacancy of 
the impurity to the Inn position. Since this process pro- 
ceeds spontaneously, it is not possible to find the corre- 
sponding vacancy jump frequency according to Eq. (4). 
According to the assumptions of the Le Claire and OSA 
models, the frequency of the exchange of the 5nn va- 
cancy with the atom further away from the impurity is 
considered to be the jump frequency of the Cr atom self- 
diffusion. In our DFT simulations, the 5nn vacancy will 
only slide into the 1nn position, so we can define: 
Wan = Wo, W3” = 0. (8) 
Based on all the assumptions made before, the diffu- 
sion coefficients of Cs and I in Cr matrix at 500 - 2000 K 
can be calculated according to OSA model and Le Claire 
model, respectively, as shown in Fig. 8. The results 
of Cr self-diffusion coefficients are shown for compari- 
son, which are in good agreement with the experimental 
data [51, 52]. It can be seen that the diffusion coefficients 
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of the impurity atoms are larger than the self-diffusion 
coefficients of Cr in the temperature interval considered. 
This can be attributed to the combined effect of the 
three factors analyzed in preceding paragraphs: weaker 
bonding effect by the surrounding Cr atoms, stronger at- 
traction effect for nearby vacancies and lower migration 
barriers compared to the Cr self-diffusion process. In the 
temperature range of Generation IV fast reactors around 
500 - 1000 K, the diffusion coefficients of Cs and I are 
3 - 7 orders of magnitude larger than the self-diffusion 
coefficients of Cr, which also indicates the strong pene- 
tration ability of Cs and I. The diffusion coefficients of 
all the three elements increase exponentially with the 
increasing temperature. Despite the different diffusion 
characteristics of Cs and I shown in Fig. 7, the diffusion 
coefficients of Cs and I are very close at the tempera- 
ture range considered. Because our calculations of Cr 
self-diffusion coefficients agree well with the experimen- 
tal results, the diffusion coefficients of Cs and I in BCC 
Cr from DFT calculations are reliable. 


The self-diffusion coefficient of BCC Cr, Der (cm?/s), 
can be recalculated using the follow equation based on 
parameters obtained in this study: 


4.0861 x 104 


Dey = 0.1323 x exp(— T ; 


(9) 


and the diffusion coefficient of Cs and I in BCC Cr, Des 
and Dr (cm?/s), can be fitted from Fig. 8 as below: 


(2.6915 + 0.0004) x 104 
T ) 
(10) 


Des = exp|(—5.8676+0.0051)— 


(3.4297 + 0.0016) x 104 


(11) 


Dy, = exp|(—1.8582 0.0158) — 
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Fig. 7. (Color online) Illustration of diffusion characteristics of Cs and I. (a) Spontaneous diffusion of Cs; (b) Non-spontaneous 
diffusion of I; (c) Migration barriers of I and Cr. Energies are relative to the corresponding initial states. 
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Fig. 8. (Color online) Temperature dependence of impurity 
diffusion coefficient in Cr matrix. Self-diffusion coefficients 
are given for comparison. Our results are denoted by lines 
and experimental results are denoted by dots. 


B. Impurity diffusion at Cr GB 


1. Atomic structure and electron distribution at the GB 


GB serves as a fast diffusion channel for the impu- 
rity atoms, and the existence of Cs and I atoms on the 
GB plane will change the GB structures, such as bond 
lengths and charge distributions. In this subsection, the 
impurity atom was placed at site 1 on the GB plane to 
analyze the changes and the differences between matrix 
and GB. As shown in Fig. 9, the presence of impurity 
atoms at GB does not lead to the increase of all bond 
lengths, which is different from that in matrix. Never- 
theless, the impurity atoms located at this substitution 
site significantly increase the bond length in the [210] 
direction, i.e., distancing the atoms on both sides of the 
GB plane. This facilitates the diffusion of the impu- 
rity atoms along the interface. Besides, Cs have a more 
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pronounced effect on increasing this bond length than I, 
which is in accord with the results in Fig. 3. 


Fig. 9. (Color online) Calculated bond lengths (A) at Cr GB 
for (a) pure GB, (b) GB with a substitutional Cs atom, (c) 
GB with a substitutional I atom. 


TABLE 3. Calculated Bader charge,volume and radius of 
atoms at Cr GB. The values of Cr are averaging values. Cr 
(Cs) denotes the Cr atom in Cs-doped Cr GB, and Cr (I) 
means that in I-doped one. 


Cs Cr(Cs) I Cr(I) 

Bader charge/e~ 8.962 6.001 7.600 5.992 
Atomic volume/A? 21.286 11.881 21.968 11.872 
Atomic radius/A 1.719 1.416 1.737 1.415 


The Bader charge, volume and radius [50] of atoms 
at GB were obtained using the Bader decomposition 
method [49] as in matrix, and the results are demon- 
strated in Table 3. As in the matrix, the Bader results of 
the 75 Cr atoms at GB were found close, so their average 
values were used for analysis. Comparing Table 1 with 
Table 3, it is found that the I atom gains electrons both 
in matrix and at GB, while the Cs atom loses electrons 
at GB, which is different from the case in matrix. For I 
substitution, the number of electrons gained at the GB 
is slightly less compared to that in the matrix, probably 
because the atomic structure at GB is not as compact as 
in the matrix and the electrons are transferred less be- 
tween the different nuclei. The atomic volumes of Cs and 
I at Cr GB are essentially the same as in Cr matrix. The 
difference is that the Bader volume of all atoms at GB is 
slightly larger than that in matrix, which is reasonable 
because the atomic arrangement at GB is sparser. This 


(b) 


Fig. 10. (Color online) Calculated charge density distribution (electron/Bohr*) in the (001) plane for (a) pure GB, and GB 


with a substitutional (b) Cs or (c) I atom at site 1. 


is also favorable to the migration of Cs and I atoms, and 
it can be expected that the impurity atoms on the GB 
plane will migrate more easily than in matrix. 


Both Cs and I have similar effects at Cr GB as in ma- 
trix, with slightly lower interactions around Cs than I, as 
shown in the distribution of charge densities in Fig. 10. 
It can be seen that the charge density around the im- 
purity atoms is significantly lower compared to that at 
pure Cr GB, which is similar as in Cr matrix (Fig. 4). 
The dilute charge density indicates a weak interaction 
between the impurity atoms with their surrounding Cr 
atoms, which favors the migration of the impurity atoms 
along the GB plane. Comparing the charge distributions 
in Fig. 10(b) & (c), it can be seen that the charge density 
around Cs is smaller than that around I, i.e. the Cs-Cr 
binding is weaker than the I-Cr one, which also agree 
with the change in the bond length in [210] direction in 
Fig. 9. Compared to the charge density distribution in 
the matrix (Fig. 4), the charge density distribution at 
GB is more dilute. As a result, it can be predicted that 
the diffusion rates of Cs and I at GB should be compa- 
rable, and the diffusion should be much easier than in 
matrix. 


2. Diffusion barrier and diffusion path along the GB 


To investigate the impurity diffusion behaviors along 
GB, impurity atoms were placed at interstitial site 0, 
substitutional sites 1 and 2 on the GB plane (as shown 
in Fig. 2(b)) respectively as initial structures of different 
migration processes, and the results of the optimization 
of these initial structures are shown in Table 4. It can be 
observed that the impurity atoms located in site 0 will 
spontaneously migrate to the non-occupied site 2, which 
is more stable than site 0. If a vacancy is introduced at 
site 1, the impurity atoms located at site 2 will sponta- 
neously migrate to site 1, which is more stable than site 


2. Both of these migration processes take place sponta- 431 
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Fig. 11. (Color online) Evaluation of energy during the Cs 
or I migration process between sites 0 and 1 at GB. Energies 
are with respect to the initial structure with impurity atom 
at site 1. 
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Fig. 12. (Color online) First-principles predicted paths of 


spontaneous impurity diffusion along GB plane: (a) impurity 
atoms initially located at site 0; (b) introduction of a vacancy 
at site 2; (c) spontaneous movement of impurity atoms to site 
2; (d) introduction of a vacancy at site 1; (e) spontaneous 
movement of impurity atoms to site 1. 


430 along the GB direction. 


The calculated migration barriers between site 0 and 


neously, so that the impurity atoms can readily diffuse 432 site 1 are shown in Fig. 11. The migration barriers from 
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TABLE 4. Atomic structures after ionic optimization by first-principles calculation. 


Initial structure 


Aligration: process (vacancy site, impurity site) 


Final position of I 
after optimization 


Final position of Cs 
after optimization 


site 2 > site 0 (site 0, site 2) site 2 site 2 
site 0 — site 2 (site 2, site 0) Between site 0 and site 2 site 2 
site 1 — site 2 (site 2, site 1) Between site 1 and site 2, closer to site 1 site 1 
site 2— site 1 (site 1, site 2) Between site 1 and site 2, closer to site 1 site 1 
site 1 — site 0 (site 0, site 1) site 1 site 1 
site 0 — site 1 (site 1, site 0) site 0 site 0 
(a) (b) (c) 455 IV. SUMMARY 
6 @ @© 6 @ @© 00o 
_ = m pP o © ase Diffusion behaviors of Cs and I, key elements of FCCI, 
0 @ 0 © o œ 457 in Cr coating were investigated using DFT approaches 
ee o » ee © D o0 o ass together with Le Claire and OSA diffusion models. Un- 
© © © © © © 459 der the combined effect of low binding effect with matrix 
o o0 0 o0 o o0 460 atoms and strong binding effect with nearby vacancies, 
ooo @e@ 0 6 oe eo 0 461 the migration barriers of Cs and I are both lower com- 
s 462 pared to that of Cr. As a result, the diffusion coefficients 
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Fig. 13. 
paths of impurity atoms along GB plane by overcoming a 


low migration barrier: (a) impurity atoms initially located at * 


site 0; (b) introduction of a vacancy at site 1; (c) impurity 
atoms diffusion to site 1 with a low energy barrier. 


site 0 to site 1 for Cs and I are obviously smaller com- 
pared to those in matrix (Table 2), indicating that GBs 
are fast channels for the impurity diffusion in Cr. Nev- 
ertheless, it can be seen that Cs or I needs to overcome a 
high energy barrier (up to 4.90 eV or 3.53 eV) to migrate 
from site 1 to site 0, meaning that impurity diffusion is 
directional in some GB regions. 


In summary, two more possible paths were predicted 
for the diffusion of impurity atoms along GB plane, as 
shown in Fig. 12 & 13. Both migration paths require the 
introduction of vacancies at the GB plane, and lots of 
vacancies exist at the GB region of real materials. Fur- 
thermore, the vacancy concentration can be significantly 
increased by the combined effect of high temperature and 
irradiation under service conditions on the inner surface 
of the cladding in nuclear reactor. In Fig. 12, vacancies 
at sites 2 and 1 were introduced in succession, and the 
impurity atoms at site 0 move spontaneously, sliding se- 
quentially to site 2 and then to site 1. In Fig. 13, one 
vacancy was introduced at site 1, and due to the low mi- 
eration barrier from sites 0 to 1, the impurity atoms can 
easily move to site 1. 


(Color online) First-principles predicted diffusion , . 
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of Cs and I are both 3-7 order of magnitude larger than 
that of Cr’s self-diffusion coefficient below 1000 K, which 
corresponds to the temperature range of Generation IV 
fast reactors. Despite the difference of migration barri- 
ers and displacements of the elemental diffusion process 
of Cs and I, and also the different models used, the dif- 
fusion coefficients obtained for Cs and I are basically 
the same order of magnitude at temperature considered. 
Meanwhile, the results of Cr self-diffusion agree well with 
experimental values, judged by which our results can be 
considered accurate. Based on our results, the diffusion 
paths of Cs and I in the inner-surface Cr coating of the 
cladding were predicted. Cs and I are more likely to dif- 
fuse along the GB, and the intergranular corrosion still 
requires for further investigations. 


V. APPENDIX 


TABLE 5. List of abbreviation. 


Abbreviation Definition 
BCC Body Centered Cubic 
cNEB climbing image Nudged Elastic Band 
DFT Density Functional Theory 
FCCI Fuel Cladding Chemical Interaction 
GB Grain Boundary 
nn nearest neighbors 
OSA Oversized Solute-Atom 
PAW Projector Augmented Wave 
PBE Perdew-Burke-Ernzerhof 
VASP Vienna ab initio simulation package 
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